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Abstract 

We examine the rare decays B s — > and B — > Kl + l~ in the framework of the 
topcolor-assisted technicolor (TC2) model. The contributions of the new particles 
predicted by this model to these rare decay processes are evaluated. We find that 
the values of their branching ratios are larger than the standard model predictions 
by one order of magnitude in wide range of the parameter space. The longitudi- 
nal polarization asymmetry of leptons in B s — ► l + l~ can approach O(10" 2 ). The 
forward-backward asymmetry of leptons in B — > Kl + l~ is not large enough to be 
measured in future experiments. We also give some discussions about the branching 
ratios and the asymmetry observables related to these rare decay processes in the 
littlest Higgs model with T-parity. 
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I. Introduction 

The study of pure leptonic and semileptonic decays of B meson is one of the most 
important tasks of B physics both theoretically and experimentally. These rare B decays 
are sensitive to new physics (NP) and their signals are useful for testing the standard 
model (SM) [1J. So far, a lot of works have been concentrated on these decays. In the 
SM, there are no flavor changing neutral current (FCNC) processes at the tree level 
and the leading contributions to these decays come from the one-loop level. So these rare 
decays are rather sensitive to the contributions from the NP models beyond the SM. 
Studying of the observables of the asymmetries, such as the CP asymmetry [2], longitudi- 
nal polarization (LP) asymmetry Alp [3], and forward-backward (FB) asymmetry Afb 
[I] etc, interests experiments in testing NP. Certainly, their detection requires excellent 
triggering and identification of leptons with low misidentification rates for hadrons. The 
precision measurement needs further studying. 

The quark level transition b — > sl + l~ is responsible for both the purely leptonic decays 
B s — > l + l~ and the semileptonic decays B — > Kl + l~(l = e,fi,r). The decay B s — > 
will be one of the most important rare B decays to be studied at the upcoming large 
hadron collider (LHC), and so far the upper bound on its branching ratio is [5] 

Br(B s -> n + iT) < 5.8 x i(T 8 (95% C.L.). (1) 

The branching ratios of B — > Kl + l~ observed by BaBar collaboration and Belle collabo- 
ration are [HI E] 

Br(B -> Kl + l~) = (5.7+H) x KT 7 , (2) 

which is close to the SM prediction [HE]- However, due to the errors in the determina- 
tion of the hadronic form factors and the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
element \Vt s \, there is about 20% uncertainty in SM prediction. The experimental mea- 
surement values of rare decay processes B s e + e~ , t + t~ will be discussed later. 

We also consider other observables of the purely leptonic and semileptonic decays for 
the B meson, which are sensitive to scalar/pesudoscalar new physics (SPNP) contribu- 
tions to b —>■ s transitions. They are forward-backward asymmetry Aps of leptons in 
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B — > Kl + U and longitudinal polarization asymmetry App of leptons in 5 S — > . The 
observable App was introduced in Reference [3], though the corresponding analysis in the 
context of K — > yU + /i~ had been carried out earlier [9]. The average App in the rare decay 
processes B — > Kl + l~ has been measured by BaBar collaboration as [6J 

(A FB ) = 0.151°;^ ±0.08. (3) 

This measured value is close to zero and has a high experimental error. As the values 
of Alp and App predicted in the SM are nearly zero, any nonzero value of one of these 
asymmetries is a signal for NP. This is the main reason we focus on these observables. 

In literature, there are numerous studies of the quark level decays b — > sl + U both in 
the SM and in some NP models. Recently, Reference [TUl [TTj have studied the sensitivity 
of these rare decay processes to the radius R in the universal extra dimension (UED) 
model. In the supersymmetry (SUSY) models, extensive works have been taken to the 
branching ratios of these rare decays, and some of these discussions are related to the 
asymmetry aspect p2J HB]- These decays have also been discussed in the littlest Higgs 
model with T-parity (called the LHT model) [Hj, they have verified that the LHT model 
can enhance the branching ratios of these decays [15J. However, they have not discussed 
the asymmetry observables, we will give some discussions on these observables in the 
framework of the LHT model. 

In the framework of the top color-assisted technicolor (TC2) model [16], Reference 
[TTj has calculated the branching ratios of quark level b — > sl + l~ decays. They consider 
the contributions of the non-universal gauge boson Z' predicted by this model. Their 
numerical results show that the enhancement is quite large when the mass of Z' is small. 
Reference [TBJ has calculated the contributions coming from the pseudoscalar top-pions 
predicted by this model to the branching ratios of the decays B s — > l + l~ . Reference [JJJ] 
has evaluated the contributions from both the neutral and charged scalars predicted by 
this model, the branching ratios can be enhanced over the SM predictions by two orders 
of magnitude in some part of parameter space. So far, we have not seen the study of 
the asymmetry observables for these two decays in the framework of the TC2 model, 
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and furthermore the former discussions on the branching ratios have not considered the 
contributions induced by all the particles predicted by this model. 

In this paper, we consider the contributions coming from all of the new particles 
predicted by the TCI model to the branching ratios and asymmetries related to the 
rare decay processes b — ► sl + F . Compared with the predictions in the SM, our results 
show that the contributions to the branching ratios and the asymmetries come from two 
aspects. First, the Wilson coefficients of these processes receive additional contributions 
from the non-universal gauge boson Z' and charged top-pions. Second, the neutral top- 
pion and top-Higgs can give contributions through newly introduced scalar /pesudoscalar 
operators. For comparison, we also give our results in the LHT model, considering 
different parametrization scenarios. 

This paper is arranged as follows. In the following section, we will summarize some 
elementary features of the TC2 model. In Sec. Ill we present our calculation on the decay 
processes B s — > l + F . The decay processes B — > Kl + F will be studied in the Sec. IV. In 
Sec. V we give simple discussions on the above questions in the LHT model. Conclusions 
are given in Sec. VI. 



The TC2 model [16J is one kind of the phenomenological viable models, which has all 
essential features of the topcolor scenario. The TC2 model generates the large quark mass 
through the formation of a dynamical ti condensation and provides possible dynamical 
mechanism for electroweak symmetry breaking (EWSB). The physical top-pions(7r° ,=t ), 
the non-universal gauge boson (Z') and the top-Higgs (h®) are predicted. The presence 
of the physical top-pions ir®^ in the low energy spectrum is an inevitable feature of the 
topcolor scenario, regardless of the dynamics responsible for EWSB and other quark 
mass. The flavor-diagonal (FD) couplings of top-pions to fermions can be written as 



II. The TC2 model 
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where m* t = m t (l — e), v w = v/y2 = 174 GeV, F n « 50 GeV is the top-pion decay 
constant. The ETC interactions give rise to the masses of the ordinary fermions including 
a very small portion of the top quark mass, namely emt with a model dependent parameter 
£ < 1, and = m& — 0.1em t [23- The factor — - reflects mixing effect between 



top-pions and the Goldstone bosons. 

For the TC2 model, the underlying interactions, topcolor interactions, are non-universal 
and therefore do not posses Glashow-Iliopoulos-Maiani (GIM) mechanism [26J. One of 
the most interesting features of ir®^ is that they have large Yukawa couplings to the third- 
generation quarks and can induce the tree- level flavor changing (FC) couplings [271 128] . 
When one writes the non-universal interactions in the quark mass eigen-basis, it can in- 
duce the tree- level FC couplings. The FC couplings of top-pions to quarks can be written 
as [TTJETJ: 

%L ^ ~ - \iK% R K«* L t L c R 4 + y^Kh c * R K% L c R b L 7rt + ^2K^ R K^ L b L c R ^ 

V 2t V V w L 

+y/2K^ R KS L t R a L i r f + ^2K^ R K s £' L s L t R ^\ , (5) 



where Kul(r) an d K^uja are rotation matrices that diagonalize the up-quark and down- 
quark mass matrices M v and M Dl i.e., K^ L M V K UR = M$ ia and K+ L M D K DR = M% a , 
for which the CKM matrix is defined as V = K^ l Kdl- To yield a realistic form of the 
CKM matrix V, it has been shown that the values of the coupling parameters can be 
taken as 



Kl « KS l « Kg L « 1, < v/27^. (6) 



In the following calculation, we will take K\j R = \j2e — e 2 and take e as in the range of 
0.03 — 0.1 [TO]. The TC2 model predicts the existence of the top-Higgs h®, which is a ti 
bound and analogous to the a particle in low energy QCD. It has similar Feynman rules 
as the SM Higgs boson, so we don't list them. 

Another significant feature of the TC2 model is the existence of non-universal gauge 
boson Z', which may provide significant contributions to some FCNC processes because 
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of its FC couplings to fermions. The FC b — s coupling to Z' can be written as 

Cz? = -y cot 9'Z>» | l -D bb D b £*s Ll ,b L - 2 -D bb D b £*s Rl ,b R + h.c. j , (7) 

Dl,Dr are matrices which rotate the down-type left and right hand quarks from the 
quark field to mass eigen-basis. The FD couplings of Z' to fermions, which are relative 
to our calculation, can be written as [TBI El EH ED] : 



1_ 1- 1- 2_ 

-TffTL ~ TR'fTR + -ztlSftL + -^> L ^b L + -t R ^t R 

2 o o 6 



- tan 2 9'Z' 



\s L ^s L - \s R ^s R - ^lY/^l ~ UrI^Hr 
o 6 2 



-e L Ye L ~ e R ^e R 



where K\ is the coupling constant and 9' is the mixing angle with tan 6' 
the ordinary hypercharge gauge coupling constant. 



(8) 
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In the following sections, we will use the above formulae to calculate the contributions 
of the TCI model to the rare decay processes B s — > and B — > Kl + l~ . 
III. The contributions of the TC2 model to the rare decay pro- 
cesses B s — > l + l~ 

The TC2 model can give contributions to rare B decays two different ways, either 
through the new contributions to the Wilson coefficients or through the new scalar or 
pseudoscalar operators. The most general model independent form of the effective Hamil- 
ton for the decays B s — > including the contributions of NP has the form: 



H(B, -> l + r) =H + H 1 



with 



(V t * s V tb ) <R A (s W b)(lj^ 5 l) 



aG} 
V2n' 



{V tb V t * s )\R s (sP R b) (11) + R P (sP R b) (h 5 l) 



(9) 



(10) 



(11) 



Where Hq represents the SM operators and Hi represents the SPNP operators. Here 
Pl,r = (1 T 75)/2, Rs, Rp and Ra denote the strengths of the scalar, pseudoscalar, and 
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axial vector operators, respectively [21]. In our analysis we assume that there are no 
additional CP phases apart from the single CKM phase, thus Rs and Rp are real. In 
the SM, the scalar and pseudoscalar couplings Rs and Rp receive contributions from the 
penguin diagrams with physical and unphysical neutral scalar exchange and are highly 
suppressed to 0(lCr 5 ). The coupling constant of the axial vector operator Ra can be 
expressed as Ra = Y SM (x)/sm 2 9 W , where Y SM (x) is the SM Inami-Lim function [32] . 
which has been listed in Appendix A. These coupling constants will receive contributions 
coming from the non-universal gauge boson Z' and the scalars rf^, h®. 
A. The contributions of the nonuniversal gauge boson Z' 



b b t s b b t s 




b t s s b t s s 





Figure 1: Penguin diagrams of Z' contributing to B s — ► / + / in the TCI model. 

In the TC2 model, the non-universal gauge boson Z' can give corrections to the SM 
function Y(x), which directly determine the coupling constant Ra- The relevant Feyn- 
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man diagrams have been shown in Fig{TJ In these diagrams, the Goldstone boson <fi is 
introduced by the 't Hooft-Feynman gauge, which can cancel the divergence in self-energy 
diagrams. Because the couplings of Z'WW, Z'cjxj) and Z'W<p do not exist in the TC2 
model, the diagrams that including the above couplings are not present. The small inter- 
ference effects between Z' and Z are not considered here. In this situation, the function 
Y TC (xt) for / = e, /i is obtained as follows: 

Y TC (xt) = l^^(C ab (x t ) + C c (x t ) + C d (x t )), (12) 

M z , 

here Xt = m* t 2 /M^y- The factor —tan 2 6' does not exist for the decay process B s — > t + t~ 
which can be seen from Eq. (El). The formations of C a b(x t ), C c (x t ) and Cd(x t ) can be 
easily obtained in the framework of the TCI model using the method in Reference [22] • 
The detailed expression forms of these functions are listed in the Appendix B. 




Figure 2: Tree level diagram of Z' contributing to B s — > l + l within the TCI model. 

The non-universal gauge boson Z' has FC coupling with fermions as shown in Eq. (J7J) , 
the tree level Feynman diagram contributing to the decay processes B s — > l + l~ has been 
shown in Figj2j The contributions can be obtained by directly calculating Figj2] using the 
standard method in Reference [29J, and the B s width can be written as: 

r < B - - i+r » - j^A^f^^ 1 ^ 2 ^ (it) 1 ' (13) 

where 

5 bs = D b L b D b L s * + 2D%D b £*. (14) 

X(6') = cot 8' for / = r, and X(9') = tan^' for I = e and /i. fs s is the decay constant of 
B s meson. 
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B. The contributions of the scalars (71"°' , 



b b t s b b t s 




Figure 3: Scalar particles contributing to B s — ► / + / in the TC2 model. 



The scalars predicted by the TC2 model give contributions to the decay processes 
B s — > Z + /~ through corrections to the coupling constants in Eq. (TTOl) and Eq. ( TTTT) . The 
relevant Feynman diagrams are displayed in FigSJ, in which (a) shows the contributions 
of neutral top-Higgs h® and top-pion 7r° to the couplings Rs and Rp, respectively; (6), 
(c) and (d) show the contributions of the charged top-pions 717" to the coupling Ra- The 
expression of the coefficient Rs can be written as: 

= V^ER ( <W C (x t ) + VtsmirntmfMjy \ 

\2V2sin 2 9 w F^ml 0t ^vg\F^m\, K sJ J 

Here x s = m* t 2 /M|, Ms is the mass of the top-pions and g% is the SU(2) coupling constant. 
C(xt) is the Inami-Lim function in the SM [32J. Since the neutral top-Higgs coupling with 
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fermions is different from that of neutral top-pion by only a factor of 75, the expression 
of Rp is same as that of Rs except only for the masses of the scalar particles. In our 
numerical estimation, we will take m^o = m h o = Ms- In this case, Rp = Rs- 

The charged top-pions 7r =l= give contributions to the SM function Y(x) via the diagrams 
(6), (c) and (d) in FigJSl the expression of the function Y TC '(x s ) can be written as: 



Y 
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(16) 



C. Numerical results 

The branching ratios of the decay processes B s — ► can be written as[3]: 



Br{B s ^l + r) 
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64tt 3 

Here tb 3 is the lifetime of B s . 

The longitudinal polarization asymmetry of the final leptons in B s — > / + /~ is defined 
as follows [3]: 



4 ± [T(g f -, g;+) + r(=Fg f -, ±s f +)] - [r(±s t -, T-s;+) + r(-s;-, -si+)] , , 

iP " [r( S/ -, Si+ ) + r(^-,± S/+ )] + [r(± S/ -,^ + ) + r(- S/ -,- s , + )] ' 1 j 

s;± are defined into one direction in dilepton rest frame as (0, i-r^r). For only one direc- 
tion, there are no difference between the final leptons, thus there is A\ p = A~[p = Alp. 
Then the Alp can be written as: 



a lp {b s -> i+r 
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(20) 



(B s -> /+/-) ~ because i2 s ~ O(10" 5 ) in the SM. 
Before giving numerical results, we need to specify the relevant SM parameters. These 
parameters have mainly been shown in Tabled! We take the coupling constant K±, the 
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= 1.166 x 1(T 5 GeV~ 2 


™b s ~- 


= 5.366 GeV 


a = 


7.297 x 1(T 3 


m B = 


5.279 GeV 


r B s - 


= (l-437iffi) x 10" 12 . 


v tb = 


1.0 


r B d 


= 1.53 x 10~ 12 s 


v ts = 


(40.6 ±2.7) x 10~ 3 


«V ■ 


= 0.105 GeV 


IBs = 


(0.259 ±0.027) GeV [33] 


M w 


= 80.425(38) GeV 


sin 2 e w = 0.23120(15) 



Table 1: Numerical inputs used in our analysis. Unless explicitly specified, they are 
taken from the Particle Data Group |34j . 

model dependent parameter e, the mass of non-universal gauge boson Mz> and the mass 
of scalars Mg as free parameters in our numerical estimation. The value of Mg remains 
subject to large uncertainty [20] • However, it has been shown that its value is generally 
allowed to be in the range of a few hundred GeV depending on the models [21]. In our 
numerical estimation, we will assume that Mg is in the range of 200GeV ~ 500GeV. 
The lower bounds on Mz> can be obtained from dijet and dilepton production in the 
Tevatron experiments [22] or BB mixing |23| . However, these bounds are significantly 
weaker than those from the precision electroweak data. Reference [21] has shown that, to 
fit the precision electroweak data, the Z' mass Mz> must be larger than 1 TeV. In our 
numerical estimation, we will assume that the values of the free parameters e, Ki and 
M z < are in the range of 0.03 ~ 0.1, ~ 1 and 1000 GeV ~ 2000 GeV, respectively. 

First we give our numerical results of the decay processes B s — ► l + l~ induced by the 
non-universal gauge boson Z'. The branching ratios of B s — > l + l~ are plotted in FigJUas 
function of the mass parameter Mz> for K\ = 0.4 and 0.8, in which we have multiplied the 
factors 10 7 and 10 3 to the values of Br(B s — > e + e~) and Br(B s /i + /i~), respectively. 
From these figures one can see that the values of Br(B s — > r + r~) are sensitive to the 
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Figure 4: The branching ratios of B s — > l + l as function of the parameter Mz> 
for Kx = 0.4 (a) and K x = 0.8 (b). 
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mass of Z', they increase as the mass parameter Mz> decreasing. For I = e,fi, the 
values of their branching ratios are not so sensitive to the parameter Mz>- Because the 
contributions of Z' to Br(B s — > e + e~) and Br(B s — > fi + fi~) are small relative to the SM 
contributions. The values of the corresponding branching ratios are both below 0(1O -9 ) 
which are not easy to be observed in current collider experiments. The contributions of 
Z 1 to the branching ratio of the decay B s — > t + t~ are large, since the non-universal gauge 
boson Z' has large couplings to the third generation fermion with respect to the first two 
generations, it can make the branching ratio value reach 0(1O -6 ) with reasonable values 
of the free parameters. 

The branching ratios of B s — > l + l~ contributed by the scalars (7r° ,=l: and h®) are plotted 
in FigJ3] as function of the mass parameter Ms for e = 0.04 and 0.08, in which we have 
multiplied the factors 10 7 and 10 2 to the values of Br(B s — > e + e~) and Br(B s fi + fi~), 
respectively. It is obvious that the values of the branching ratios for these decays increase 
as the parameter Ms decreasing. Furthermore, the enhancement to the branching ratio 
of the decay process B s —>■ is larger than that of the Z' contributions by an order 

of magnitude. 
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Figure 5: The branching ratios of B s — > l + l 
e = 0.04 (a) and e = 0.08 (b). 
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Figure 6: The longitudinal polarization asymmetry in B s — > Z + / as function of the 
parameter M s for e = 0.04 (a) and e = 0.08 (b). 



13 



The value of Br(B s — > e + e~) is smaller than that of Br(B s — > by five orders 

of magnitude, which is because it is suppressed by ml/m 2 ^ with respect to fi channel. 
The branching ratio for t + t~ mode is enhanced by a factor of 10 2 to fi channel, its value 
can reach O(10~ 6 ) by our calculation. However, the r + r~ channel is still not easy to be 
observed under present experimental precision, while the current experimental upper limit 
for Br(B s -> r+r - ) from the BARBAR collaboration is 4.1 x 1(T 3 at 90% C.L. [35]. So the 
experimental searches for B s — > have focused on the \i channel, and we only discuss 
this channel. Comparing with the SM prediction Br(B s — > / u + ^) = 3.86 ± 0.15 x 10~ 9 
PP, the contributions of the new scalars predicted by the TC2 model can enhance this 
value by one order of magnitude, so our results are more approach to the experimental 
data given by Eq. (PQ). 

Obviously, the non-universal gauge boson Z 1 has no contributions to the SPNP opera- 
tors, so it was not considered in this subsection. The longitudinal polarization asymmetry 
Alp contributed by the new scalars predicted by the TC2 model as function of the param- 
eter Ms are plotted in FigJHl From these figures one can see that the Alp is sensitive to 
the mass of the scalars, especially for I = fi,r, however it is less sensitive to the parameter 
e. The values of the asymmetry Alp can reach nearly 4% for I — fi,r when the mass of 
the scalars get to 200 GeV. 

IV. The contributions of the TC2 model to the rare decay pro- 
cesses B — > Kl + l~ 

The effective Hamilton for the decay B — > Kl + l~ is similar to that of B s — > as 
shown in Eq. which is constituted by two parts. The SPNP part is same as the 
expression shown in Eq. ( TlTl) . In the framework of the TC2 model, The Hq part can be 
written as I3T1: 



H = ^V tb V: s \cf{- S1 ^PLb) h,l + C 10 {s^P L b) h^ld 
- 2— ^-m b (sia^Ppb) l^lj . 



(21) 



Here g M is the sum of 4-momenta of Z + and / . The Wilson coefficients , Cg^ and 
Cio contain two parts of contributions from the SM and the TC2 model. 
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Similar to the decay processes B s — > Z + /~, the non- universal gauge boson Z' give 
contributions to the Wilson coefficients Cg and Cio, the relevant Feynman diagrams are 
same as Fig {U and the relevant functions Y TC (xt) and Z TC (a;t) have same expressions as 
shown in Eq. (1X2*1) . 

The charged top-pions 7r^ can give contributions to the Wilson coefficients C% and 
Cg . The relevant Feynman diagrams are similar to FigJHl However, these penguin 
diagrams are induced by 7 penguins, gluon penguins and chromomagnetic penguins. The 
coefficients Cy and Cg^ can be expressed in terms of the corresponding functions Di(x s ), 
Ei(x s ) and E[(x s ), which are added to the corresponding SM functions D (x t ), E (x t ) 
and E' (xt) [32] • The detailed expression forms of the these functions are [36] : 



We can obtain the corrected Wilson coefficients Cj , Cg and C w with these corrected 
functions using the relevant expressions of these coefficients in References [TQl [SH] , which 
are listed in Appendix C. The neutral top-pion 7if and top-Higgs can also give contri- 
butions to these decay processes through the SPNP operators, and the expression forms 
of Rs {Rp) are same as those shown in Eq. (fToT) . 

The branching ratios Br(B — > Kl + l~) (I = e,fi and r) contributed by the gauge 
boson Z' are plotted in FigfT] as a function of the mass parameter Mz> for two values 
of Ki, in which we have multiplied the factor 1CT 1 and 10~ 2 to the branching ratios of 
decays B — > Kfi + fi~ and B —>■ Kt + t~ respectively. From this figure one can see that the 
values of the branching ratios for I = e, fi and r increase as the parameter Mz> decreasing. 
However, the branching ratios for / = e are not sensitive to the parameter Mz' as shown 
in these figures. The values of the branching ratios for / = e and fi are not sensitive 
to the parameter K x . For K x = 0.4 and lOOOGeV < M Z / < 2000GeV, the values of 
Br(B — > Ke + e~) and Br(B — * K{i + [i~) are in the range of 6.1 x 1CT 8 ~ 4.4 x 1CT 8 and 




(22) 



(23) 



(24) 
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Figure 7: The branching ratios of B — > Kl + l~ as function of the parameter Mz> 
for K x = 0.4 (a) and K x = 0.8 (b). 

3.0 x 10~ 7 ~ 1.2 x 10~ 7 , respectively. 

The branching ratios of the decay processes B — > Kl + l~ contributed by the scalars 
(7r°' , h®) are plotted in FigJH] as function of the mass parameter Ms for e = 0.04 and 
0.08, in which we have multiplied the factors 10 _1 to the branching ratio of B — > Kfi + fi~ . 
From these figures, one can see that the values of the branching ratios of these decay 
processes increase as the parameter Ms decreasing. All of their values are not sensitive to 
the parameter e. The contributions of the scalars for I = e and [i are comparable to those 
of the non-universal gauge boson Z\ the values of the branching ratios of B — * Ke + e~ 
and B — > Kfi + jjT contributed by both the scalars and the non-universal gauge boson 
can reach (9(10 -7 ), which give an explanation to the deviation between the experimental 
data and the SM predictions in Reference [8]. While the scalar's contribution to the 
decay process B — > Kt + t~ is smaller than that of the non-universal gauge boson Z' 
by two order of magnitude and therefore can be neglected. When the Z' mass is in the 
range of 1000 GeV ~ 2000 GeV, the values of BriB — > Kt + t~) are in the range of 
7.0 x 10~ 6 ~ 1.7 x 10~ 6 . This result is 2 orders of magnitude larger than the e and \i 
channel, which is because of the large coupling of Z' to the third generation fermions. 
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M s (GeV) 



(a) e = 0.04 



(b) e = 0.08 



Figure 8: The branching ratios of B — > Kl + l as function of the parameter Ms for 
e = 0.04 (a) and e = 0.08 (b). 



500 



The normalized forward-backward (FB) asymmetry can be defined as [3]: 

JO aCUbU dzdcosB J-l UUJbu t 



A FB (z) 



dzdcosd 
~~cPT 



Jo dc0s9 diS^ + /-I dc0s9 dzdcos8 

After the integral calculation of FB asymmetry gives, 

2T B Y Q m l l3 2 R s J dza x {z)4){\,k 2 ,z) 



(25) 



(A FB ) 



Br(B -> fa+/- 



(26) 



where is the lifetime of B meson and Br(B — > Kl + l ) is the total branching ratio of 
B — > Kl + l~ and Tq is the total width of the S meson, which can be written as: 



r 2 ™ 2 



ai (z) = ~(1 - fc 2 )C 9 /oW/ + (*) + (1 - k)C 7 fo{z)f T (z) . 



(27) 



(2? 



Other relevant functions such as 0(1, fc 2 , 2) are listed in Appendix C. The form factors 
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Figure 9: In the TCI model, the forward-backward asymmetry in B — > Kl + l as 
function of Ms for the parameter e = 0.04 (a) and e = 0.08 (b). 



/+, fo and J't are defined in the relevant matrix elements as: 

l-k 2 . 



(K(p')\s 1 ,b\B(p)) 
(K(p')\sia tlv q v b\B(p)) 
(K(p')\sb\B(p)) 



(2p-q) fl f + (z) + (- 



(2p - q)^q 2 



(m 2 B - m 2 K ) qfl 



<U/o(*) -/+(*)] > (29) 



vtiB + m K 



msil — k 2 
rh b 



U(z) . 



(30) 
(31) 



Here, k = m K /m B , z = q 2 /m 2 3 and m b = m h /m B . The form factors /+, / an d /r can be 
calculated by using the light cone QCD approach. Their particular forms can be found 
in Reference [31J. In this paper, we assume rhb = 1- 

The production of the FB asymmetries are only sensitive to SPNP operators. From 
Eq. (I2oT) . one can see that the non- universal gauge boson Z' has no contribution to the 



FB asymmetry, so we only discuss the contributions coming from the scalars (ir t ' , h®). 

The FB asymmetry Afb of leptons in the decay processes B — > Kl + l~ are plotted in 
Figj9]as function of the parameter Ms for e = 0.04 and 0.08, in which we have multiplied 
the factors 10 5 and 10 to the value of Ap B (B — > Ke + e~) and Ap B (B — > K^i + ^r) respec- 
tively. From this figure one can see that the value of Afb is smaller than O(10~ 3 ) in most 
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of the parameter spaces. Comparing its experimental measurement value, this value is 
not large enough to be observed in experiments. One can see that the contributions of 
the TC2 model to the FB asymmetry in these decay processes are smaller than those of 
the SUSY models. Considering the uncertainty in measurements, it is very difficult to 
detect the signals of the TC2 model through measuring the FB asymmetry about these 
decay processes. 

V. The contributions of the LHT model to the rare decay pro- 
cesses b —> sl + l~ 

The LHT model [H] is based on an SU(5)/SO(5) global symmetry breaking pattern. 
A subgroup [ST/ (2) x U(l))i x [SU(2) x U(l)] 2 of the SU(5) global symmetry is gauged, 
and at the scale / it is broken into the SM electroweak symmetry SU(2)l x U(l)y 
T-parity is an automorphism which exchanges the [SU(2) x C/(l)]i and [SU(2) x £/(l)]2 
gauge symmetries. The T-even combinations of the gauge fields are the SM electroweak 
gauge bosons and B^. The T-odd combinations are T-parity partners of the SM 
electroweak gauge bosons. 

After taking into account EWSB, at the order of i> 2 /f 2 , the masses of the T-odd set 
of the SU(2) x 17(1) gauge bosons are given as: 

Mb h = ^[1-^}, M ZH ^M WH = gf[l-^], (32) 

where / is the scale parameter of the gauge symmetry breaking of the LHT model, g' is 
the SM U{l)y gauge coupling constants. Because of the smallness of g', the T-odd gauge 
boson Bh is the lightest T-odd particle, which can be seen as an attractive dark matter 
candidate [37] . To avoid severe constraints and simultaneously implement T-parity, it is 
necessary to double the SM fermion doublet spectrum [14, 38J. The T-even combination is 
associated with the SU(2)l doublet, while the T-odd combination is its T-parity partner. 
The masses of the T-odd fermions can be written in a unified manner as: 

M Fi = V2hf, (33) 

where fc, are the eigenvalues of the mass matrix k and their values are generally dependent 
on the fermion species i. 
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The mirror fermions (T-odd quarks and T-odd leptons) have new flavor violating 
interactions with the SM fermions mediated by the new gauge bosons (Bjj, W^, or Zh), 
which are parametrized by four CKM-like unitary mixing matrices, two for mirror quarks 
and two for mirror leptons 



V Htt , V Hd} V Hh V Hv . (34) 

They satisfy: 

Vff u VHd = Vckm, V^ u Vhi = Vpmns- (35) 

Where the CKM matrix Vckm is defined through flavor mixing in the down-type quark 
sector, while the PMNS matrix Vpmns is defined through neutrino mixing. 

The contributions of the LHT model to the rare decay processes b — > sl + l~ are mainly 
coming from the corrections to the Wilson coefficients, which related to the SM Inami- 
Lim functions [32] . The branching ratios of the decay processes B s — > in the SM 
depend on a function Y$m and the LHT effects enter through the modification of the 
function Ysm [HH]. With the LHT effects Ysm is replaced by [T5] : 



Y S = Y SM + Y e ™ + ^_, (36) 

where "p even and Y° dd represent the effects from T-even and T-odd particles, respectively. 
The branching ratios normalized to the SM predictions are then given by: 



Br(B s -> m- 



Y, 



SM 



2 



(37) 



Br (B s -> Z+/-) SM 

which Br(B s — > / + Z~)sm are the branching ratios predicted by the SM. Their particular 

numerical values of the branching ratios for the decay processes B s — ► in the LHT 
model are listed as follows: 

Br(B s ^e + e-) = (1.36 ± 0.05) x 10" 13 , (38) 

Br{B s ^n + n~) = (5.79 ± 0.23) x 10" 9 , (39) 

Br(B s — ► r + r _ ) = (1.23 ± 0.05) x 10" 6 . (40) 



20 



The branching ratios of the decay processes B — > Kl + l~ in the 5M depend on the 
functions Y$m, Zsm and D' (x t ) (D' (x t ) is same as in B — ► X s 7 [15]), the LifT effects 
enter through the modification of these functions. The modifications of the function Ysm 
has been given above, and the modifications of the function Zsm is given by [T5"j I3"9"] : 

Z s = Z SM + ^ even + -f^r, (41) 

Ar 

where Z even and represent the effects coming from T-even and T-odd particles, re- 
spectively Similar with Sec. IV, we can calculate the contributions of the LHT model 
to the decay processes B — > Kl + l~ . With reasonable values of the free parameters in the 
framework of the LHT model, the maximum values of the branching ratios for the rare 
decays B —>■ Kl + l~ are: 

Br(B -> Ke + e~) = 9.66 x 1(T 6 , (42) 
Br(B -> Kfi+fi-) = 6.56 x 10" 6 , (43) 
Br(B -> Kt + t-) = 2.99 x 10" 7 . (44) 

These numerical results are obtained by calculating the relative correction to the SM 
predictions in the framework of the LHT model, while the SM predictions exist the 
uncertainty coming from the next-to-leading logarithmic (NLO) contributions and the 
long-distance contributions, for which the Br(B — > Kl + l~) are a little disparity away 
from their respective experimental upper limits [41] . However, there is no disagreement 
with experiment in some parameter ranges while the corrected effects is no more than 15 
percent. 

The contributions of the LHT model to the asymmetry observables Aps and A^p in 
the rare decay processes b —± sl + l~ mainly come from the new neutral scalar particles. 
For the B s meson, there is an unitarity relation of the VhcL matrix |39j : 

d s) + # + d S) = 0, (45) 

where Q = V^V^ d . Considering this relation, the calculations of the relevant Feynman 
diagrams similar to Figj3] equal to zero. Hence, in the framework of the LHT model, the 
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total contributions induced by the neutral scalars equal to zero. The contributions to the 
A FB and A L p is close to the predictions in the SM. 
VI. Conclusions 

The SM is a very successful theory but it can only be an effective theory below some 
high energy scales. To completely avoid the problems arising from the elementary Higgs 
field in the SM, various kinds of dynamical electroweak symmetry breaking models have 
been proposed, among which the topcolor scenario is attractive because it can explain the 
large top quark mass and provide a possible EWSB mechanism. The TC2 model has all 
essential features of the topcolor scenario. It is expected that the possible signals of the 
TC2 model should be detected in the future high energy collider experiments. 

In this paper we consider the contributions of the TC2 model to observables related 
to the decay processes B s — > l + l~ and B — > Kl + l~ . We find that the TC2 model can 
enhance the branching ratios of the SM predictions for these decay processes B s — > / + /~ 
and B — > Kl + l~ . In wide ranges of the free parameter space, it is possible to enhance the 
values of Br(B s — > and Br(B — > Kl + l~) by one order of magnitude. In the TC2 

model, the non-universal gauge boson Z' gives main contributions to Br(B s — > t + t~), 
while the contributions of Z' to Br(B s — > e + e~) and Br(B s — > are comparable with 

those of the new scalars (71^, h° t ). For the decay processes B — > Ke + e~ and B — > Kfi + /i~ , 
the contributions of Z 1 are comparable with those of the scalars. While the contributions 
of the TC2 model to Br(B — > Kt + t~) mainly come from Z' . 

The production of the asymmetries are only sensitive to SPNP operators, so there are 
no contributions of Z 1 to the relevant observables. We further calculate the contributions 
of the new scalars predicted by the TC2 model to the asymmetry observables Afb and 
Alp of leptons in the decay processes B s — > l + l~ and B — > Kl + l~ . Our numerical results 
show that, when the mass of the scalars gets to 200GeV, the values of the asymmetry 
A LP in the decay processes B s — > /x + and B s — > t + t~ can reach 4% . We hope that 
the values of Alp for I = fi,r can approach the detectability threshold of the near future 
experiments. However, the contributions of these new scalars to App are around (9(10 -4 ) 
in most of the parameter space, which are not large enough to be detected. 
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The LHT model is one of the attractive little Higgs models, which satisfies the elec- 
troweak precision data in most of the parameter space. This model can produce rich 
phenomenology at present and in future high energy experiments. New particles pre- 
dicted by this model give contributions to the branching ratios of the rare decay processes 
B s — > and B — > Kl + l~ . Reference [IS] has shown that, comparing with their SM 
predictions, the branching ratios of the decay processes B s — > l + l~ and B — > Kl + l~ can 
be enhanced by at most 50% and 15%, respectively. For comparison, we give a brief 
description and particular numerical results about these rare decays. In addition, we 
show that the neutral scalars predicted by this model can not give contributions to the 
asymmetry observables Afb and Alp- 

In conclusion, the effects of the TCI model on the branching ratios and asymmetry 
observables related to the rare decay processes b — > sl + l~ can give positive contributions 
to the SM predictions. The numerical results show that the branching ratios for these 
decays are much close to the experimental data, such as Br{B s — > The value of 

BriB — > Kt + t~) is larger than the SM prediction by one order of magnitude, which is 
hoped to be observed in the future high accuracy experiments, or the future experimental 
results may give constraints on the free parameters of the TCI model. Hence, it is 
indicated that the possible signals of the TC2 model may be observed through the above 
decay processes in future experiments. 
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Appendix 

A. Relevant functions in the SM 
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In this Appendix we list the functions in the SM that entered the present study of 
rare B decays. 



Y SM (x) 



z - 4 
x — 1 



3x 



{x-iy 



Inx 



Z SM (x t ) = --lnx t + 



9'""" ' lU(xt-l) 3 
32a; 4 - 38x 3 t - 15x 2 + I8x t 



72(x t - l) 4 



hax t , 



(46) 



(47) 



D (y) 
E (y) 

DM 

K(y) 



4 t -19y 3 + 25y 2 , y 2 (5y 2 - 2y - 6) , 
9 lny+ 36(,-l)3 + 18(,-1) 4 ln2/ ' 
2 y 2 (15-16y + V) , ^(18-lly-y 2 ) 

S lny+ 6(,-l) 4 ^ + 12(1 -y)3 
(3y 3 -2t/ 2 ) ln ^ + (% 3 + 5y 2 -7y) 



2(y-l) 4 12(y-l) 3 

3y 2 . fy 3 - 5t/ 2 - 2y) 
lny + 



(48) 
(49) 
(50) 
(51) 



2(3/ -l)* * 4(2/- 1)3 

B. Relevant functions in the TC2 model 

In this Appendix we list the functions that entered the present study of rare B decays 
in the framework of the TC2 model. 



C a b{x) 

C c (x) 
C d (x) 
C(x) 



3gi(v d + a d ) 
2/ 2 c 2 , ( 2 /:,(. r 



9l 



+ 



FJx) 



2f 2 c 



2„2 
w 



9l 



u ®u ) J 

2F 4 {x) F 5 (x) 



3(v u + a u ) 6(v 
F^x) 



(52) 
(53) 
(54) 
(55) 



-(0.5(Q-l)s 2 +0.25)' 
Here the variables are defined as: g = y/AirKi, v u>d = h~ 2Qu,dSw, s w = sin0 w , a ud = I 3 , 
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where u, d represent the up and down type quarks, respectively. 



F^x) = -(0.5(Q - 1)4 + 0.25)(x 2 ln(x)/(x - l) 2 - x/(x - 1) - x(0.5(-0.5772 

+ /n(4vr) -ln(M^)) + 0.75 - 0.5(x 2 ln(x)/(x - l) 2 - l/(x- 1)))), (56) 

F 2 (x) = (0. 5Q s 2 w - 0. 25) (x 2 ln(x) / (x - I) 2 -2xln(x)/(x -I) 2 + x/(x- I)), (57) 

F 3 (x) = -Qs 2 w (x/(x-l)-xln(x)/(x-l) 2 ), (58) 

F 4 (x) = 0.25(44/3- l)(x 2 ln(x) / (x - l) 2 -x-x/(x- 1)), (59) 



F 5 (x) = -0.25gs 2 ^(-0-5772 + /n(4vr)-/n(M^) + l-a;/n(a;)/(x-l)) 
- s 2 /6(x 2 Zn(x)/(x - l) 2 - x - x/(x - 1)). 



(60) 



C. Relevant expressions in our calculation 

In this Appendix we list the functions that entered the present study of rare B decays 
and some expressions of the relevant coefficients. 

M(B^Ki+n = f^ws 

x (K(p') \s^b\ B(p)) {Cfu(p + )^v(P-) + C w u(p + ) W v(p-)} 
-2-^m b {K{p') \sia, v q»b\ B(p)) u(p+) llx v(p_) 



+ {K{p') \sb\ B(p)) {R s u(p + )v(p_) + Rpu(P + hnv(p~)} 



(61) 



d 2 T 



dzdcosO 



G 2 F a 2 
2% 5 



|V*l£| 2 ro|0 1/a (l,* 2 , 



x 



[\A\ 2 Pl + lB^z + ^k 2 ^)^ 2 +\D\ 2 )(1- Plcos 2 9) 



+2m,(l - k 2 + z)Re(BC*) + 4m 2 \C\' 



+2rhi<p*{l,k 2 ,z)(3 ll Re(AD*) cosd 



(62) 
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A = ±(l-k 2 )f (z)R s , 

B = -7n l C 10 ^f4z)~ 1 -^(f (z)-f + (z))^ + ^(l-k 2 )f (z)R P 

C = C w f+(z), 

D = C e 9 ff f + (z) + 2C 7 



1 + k 

0(1, k 2 , z) = 1 + k 4 + z 2 - 2{k 2 + k 2 z + z) , 

Ami 2 . , . 

= -)■ (63) 

In place of CV, one defines an effective coefficient C^ e ^ which is renormalization 
scheme independent [12] : 

cf )e/ v) = #cfW) + jjfo* - ^)cf(^) + cf(^) £ (64) 

i=i 

where ?7 = , and 

CfW) = 1, CfW) = -^'(^), ^°W) = -\&{x t )\ (65) 

the superscript (0) stays for leading logarithm approximation, which is not displayed in 

the text. Furthermore: 

14 16 6 12 

ai = 23 " 2 = 23 " 3= 23 tt4 = ~23 
a 5 = 0.4086 a 6 = -0.4230 a 7 = -0.8994 a 8 = -0.1456 

h x = 2.996 h 2 = -1.0880 h 3 = -I h A = -— 

7 14 

h 5 = -0.649 he = -0.0380 /i 7 = -0.0185 h s = -0.0057. (66) 
In the Naive dimensional regularization (NDR) scheme one has 



Y(x t ) 



CM = P Q " nH + - AZ(x t ) + P E E(x t ) (67) 

S w 

where Pq DR = 2.60 ± 0.25 [12] and the last term is numerically negligible. 
Cio is fi independent and is given by 

C W = (68) 

The normalization scale is fixed to fi = fib — 5 GeV. 
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